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The impact dynamics between wet surfaces, which dominates the mechanical properties of wet 
granular matter, is studied both experimentally and theoretically. It is shown that the hysteretic 
formation and rupture of liquid capillary bridges between adjacent grains accounts reasonably well 
for most relevant cases of wet granular matter. The various dissipation mechanisms are discussed 
with particular emphasis on their relevance. Variations of the rupture energy loss with the impact 
energy are quantified and discussed. 
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The interest in granular materials among the soft- 
matter community has been strongly increasing in recent 
years. It has been widely recognized that many con- 
cepts which are well established for colloidal systems and 
glasses apply as well to granular systems [TJ [H [31 IH [S] , 
and vice versa [Tj [3 [H]. Furthermore, both colloidal 
and granular matter play a certain role as models for 
other systems which are too complex to be tractable. In 
particular, granular systems are of great interest in the 
context of dynamical systems far from thermal equilib- 
rium [iniiniiTi]. 

The term 'soft matter' applies particularly well to wet 
granulates, which can be shaped to stable structures |13j . 
but yield to rather small shear stress [T31 [TBI fTZl [15] . 
This plasticity stands out against dry granulates, such as 
the sand in an hour glass, which runs through the orifice 
like a fluid [19 . The reason for this difference is that in 
the wet system, small liquid capillary bridges form be- 
tween adjacent grains, exerting an attractive force upon 
them by means of the surface tension of the liquid [20 , 21J. 
It is clear that these only form when the liquid wets the 
material the grains consist of, which is well fulfilled for 
most sands. When a wet granulate is being sheared, or 
otherwise mechanically agitated, the repeated formation 
and rupture of the many liquid objects inside gives rise to 
considerable dissipation, which is then experienced as a 
noticeable resistance to the external drive imposed on the 
material. In order to understand the mechanical prop- 
erties of wet granular matter, it is thus indispensable 
to understand the dissipation processes connected to the 
liquid capillary bridges in detail. 

As it is well known, there is some intrinsic dissipation 
also in dry granular matter, which is responsible for the 
fact that even the perfectly dry granulate in the hour- 
glass behaves distinctly different from a regular fluid. 
Some fraction of the kinetic energy of the grains is trans- 
ferred at each impact to the microscopic degrees of free- 
dom on atomic scale. The 'heat bath' represented by 
the random center-of-mass motion of the grains, which 



easily corresponds to Giga- or even Tera-Kelvins when 
converted to temperature using Boltzmann's constant, is 
thus intimately coupled to the room-temperature heat 
bath of the atoms. Some of the most striking features 
of granular motion owe to this intrinsic non-equilibrium 
character. One usually quantifles these effects by means 
of the so-called restitution coefficient, e, which is defined 
as the ratio of the momenta before and after the impact, 
Pf = e Pi. It is found that this ratio is fairly indepen- 
dent of the initial kinetic energy of the grains in a wide 
range, although for large energies it tends to be somewhat 
smaller [5]. The most important aspect of this dissipa- 
tion mechanism is that the energy lost in the impacts, 
AEdry, scales with the impact energy, E, according to 

AEdry = (1 - e^)E (1) 

where s is approximately constant. There is thus no spe- 
cific energy scale set by this process. 

On the contrary, the liquid capillary bridges which are 
present in a wet granulate do provide their own energy 
scale. This is due to their characteristic dynamics of 
bridge formation and rupture. When two wet grains ap- 
proach each other, the liquid adsorbed on their surface 
will not react until they come into contact. At this point, 
liquid is rapidly dragged to the area of contact due to the 
interfacial forces, and a capillary bridge forms. When the 
grains withdraw from each other after the impact, the 
bridge remains intact for quite some distance, exerting 
an attractive force upon the grains. This is illustrated in 
Fig.[l^ for the idealized case of spherical grains. The an- 
gle 9 is the contact angle the liquid makes with the grain 
material, and characterizes its wetting properties. For 
complete wetting, we have — 0. At a certain critical 
separation of the grain surfaces, which depends upon the 
liquid volume of the capillary bridge, the latter ruptures 
and distributes its liquid content back onto the grain sur- 
faces. If = V/R'^ is the normalized liquid volume of a 
bridge between spherical grains of radius R, and s = s/R 
is the normalized separation of the grain surfaces, rup- 
ture occurs at 

~s. = {l + l){v'/' + 0.lV'/') (2) 
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in good approximation 
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FIG. 1: (a) A liquid capillary bridge between two spherical 
'grains' of radius R. The capillary bridge can span a certain 
distance between the grains, until it pinches off at a critical 
distance which depends upon the liquid volume in the bridge, 
(b) A typical force-vs-distance curve observed a liquid bridge 
has formed. It is shown as normalized with respect to the 
contact force, 2tt^Rcos9. 



A typical force-vs-distance curve is shown schemati- 
cally in Fig. The force is normalized with respect to 
the contact force, Fq = 27r_R7Cos6', where 7 is the surface 
tension of the liquid [22]. Since there is no liquid bridge 
(and thus no force) when the grains are approaching, the 
energy lost in the entire process of formation and rup- 
ture of the liquid bridge is given by the area under the 
descending curve, 



F{s) ds 



(3) 



For spherical grains, the shape of this curve is well 
known in the quasi-static case |22]. It corresponds to 
the force exerted by a rotationally symmetric minimal 
surface spanned between the spheres, where the liquid 
volume of the bridge and the contact angle are the main 
geometric parameters. A good approximation is 
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where S — s/yV [22] . This is in fact the curve displayed 
in Fig. [l|). It terminates at Sc = Sc/V^ as given by 



eq. Q. Using this approximation, we can evaluate the 
integral and obtain 



AEcap = 0.67 FqRVV arctan (0.35 -I- 1.68Sc 



(5) 



for complete wetting {9 — 0). 5c is at least of order 
unity, such that arctan(...) w 1.1 or larger, and it never 
goes beyond 1.57. Hence a reasonable approximation is 
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Given the fact that we have neglected side effects such 
as contact angle hysteresis, this should be as good as it 
gets. 

The presence of a defined energy loss, which does not 
scale with the impact energies of the colliding grains, has 
dramatic consequences for the collective physical proper- 
ties of the system. Most prominently, it leads to phase 
transitions which occur when the granular temperature 
is comparable to the energy loss [TU I13J |231 US]. It is 
illustrative to express the fixed energy loss in terms of 
an energy-dependent restitution coefficient. We readily 
obtain 



= A 1- 
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where AEcap is constant. Obviously, e becomes zero 
when E — AE^ap- As a consequence, an energy-driven 
phase transition occurs when the granular temperature 
comes close to AEcap [H]- This is not the case in a dry 
system described by a constant restitution coefficient. 

If the impact dynamics is sufficiently slow as compared 
to the dynamics of bridge formation and rupture, we may 
assume that the dynamic force-vs-distance curve corre- 
sponds to the quasi-static case, as represented by eq. Q. 
However, when we consider dynamical processes, as they 
take place in a sheared or otherwise agitated wet gran- 
ular material, we have to discuss the influence of this 
dynamics on AEcap- 

In doing so, we will also have to consider the energy 
loss due to viscous damping in the liquid. In reasonable 
approximation, the viscous force is given by |26j 
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for spherical grains, where rj is the viscosity of the liquid 
and V is the relative velocity of the grains at impact. If 
viscous forces are dominant, the equation of motion dur- 
ing withdrawal reads Fmsc = —ms, with the grain mass 
m, and the dot indicating the derivative with respect to 
time. Direct integration leads to 
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where vq denotes the velocity directly after the impact 
and 5 is a cutoff parameter, which may be identifled 
with the roughness of the grains [12j. vq = 9r]/8pR is 
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a characteristic velocity, where p is the density of the 
grain material. If u = In ^ , the grains stick together 
in the sense that their kinetic energy is not sufficient to 
supply the viscous energy required for reaching the rup- 
ture distance, Sc, even at zero surface tension. If v is 
considerably larger, one obtains the known result [12] 
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AE^isc ~ —R'-qvin — 



(10) 



Eq. ([9| may as well be expressed in terms of an energy- 
dependent restitution coefficient as 



= 1 - 



(11) 



where Eq — y [woln(sc/(5)]^ is the impact energy below 
which sticking occurs. We see from eqs. ^ and (111 that 
there are great similarities between the capillary and the 
viscous effects, and both strongly differ from the scale- 
free energy loss encountered with the dry systems. 

In order to investigate to what extent the concepts 
discussed above apply under realistic conditions, we have 
performed an experiment particularly designed to map 
the conditions in agitated wet granular matter as closely 
as possible. As model grains, we have used spherical 
glass beads with i? = 1 mm. These fell freely in a closed 
box containing humid air onto a wet glass plate. Initial 
heights were about 2 cm, such that the impact velocities 
were on the order of a few cm/sec, which corresponds 
to typical granular temperatures in agitated granulates. 
The motion of the glass spheres was recorded with a fast 
CCD camera, the images were subsequently analyzed by 
standard image processing techniques. 

Since all arguments put forward above concerning 
forces between two spherical grains apply as well to forces 
between a sphere and a flat wall, we have studied the lat- 
ter because of its better experimental accessibihty. We 
just have to keep in mind that in formulas developed for 
two spherical grains, the radius of the sphere must be 
multiplied by two. This corresponds to the well-known 
Derjaguin approximation, R~ff — {Ri^ + ^2^)/'^- 
though this is not an accurate expression for the system 
under study [12] , it provides a reasonable approximation 
since other effects, like contact angle hysteresis, give rise 
to larger uncertainties [T^ . 

The result is shown in Fig. [2] The height of a bounc- 
ing glass sphere is plotted as a function of time (circles) . 
On the same time scale, the total energy as obtained 
from the instantaneous velocity and height is indicated 
by the black squares. While there is some scattering in 
the vicinity of the impacts, mostly due to the finite delay 
between consecutive images, the energy is observed to be 
constant with high accuracy away from the impacts, al- 
lowing for an accurate determination of the energy level 
of each bounce. 

Below the main panel, we show three consecutive 
closeup images in the ultimate temporal vicinity of an 
impact with the glass plate. The time elapsed between 
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FIG. 2: Top: An experimental run as captured by a fast 
CCD camera (491.3 frames/sec). From the trajectories, we 
can derive the total energy as a function of time, which is 
shown as the small black squares. Bottom: three consecutive 
images taken around the time of impact with the bottom glass 
plate. The hysteretic character of the bridge is clearly seen. 
The time elapsed between two consecutive images is 2.04 ms. 



the images is 2.04 ms. The hysteretic character of the 
capillary bridge is clearly visible from its absence before 
the impact (left) and its persistence afterwards (right). 
From the difference in the energy levels of consecutive 
bounces, we can deduce the energy lost in the impact 
with the glass plate. For a quantitative analysis, we have 
to consider as well the energy loss due to viscous friction 
in the air, AEair- Using Stokes' formula, it is straight- 
forward to see that 
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where 77 is the viscosity of air. This 
low Reynolds numbers. We have Re ft 



is valid 



at 



^iJ/15.6^m, where H is the height of the bounce, such 
that Re reaches values around 30 in our experiment. 
Since turbulence sets in only at much higher Re |27j . we 
can safely assume eq. (12 1 to describe our system well. 



The most convenient way of analyzing the data is to 



4 



plot the energy of each bounce as a function of the energy 
of the previous one. Taking all dissipation mechanisms 
into account, we obtain 



e^En 



AE„ 



AE„ 



(13) 



where n numbers the bounces, and the last two terms 
depend upon En- The result is displayed in Fig.|3]i. The 
full squares represent the results for En+i as obtained 
form the experiment, while the open circles have been 
corrected for AEair according to eq. ( 12 ), which is known 



without any free parameters. As one can clearly see, the 
correction is of minor importance, as as suggested already 
by the energy data from the trajectories in Fig. [2^. From 
eq. pOl, we see that the viscous energy loss in the liquid 



is of order mvQV, which is readily checked to be below 3 
nJ in our experiments. It is thus even smaller than the 
viscous dissipation in the air, and will henceforth be ne- 
glected. The experimental error of the results displayed 
in Fig. [3^ is well below the size of the symbols. 

The three data points for the largest energies lie on a 
straight line within errors. From the slope of this line, we 
obtain the restitution coefficient connected to the solid 
impact of the glass sphere with the bottom plate. The 
result is £ = 0.82 ± 0.02. Interpreting the intercept with 
the vertical axis as the capillary bridge energy, we obtain 
AEcap = 0.07 ± 0.02/iJ. However, this intercept can be 
determined in principle for each pair of consecutive en- 
ergy levels, E„. The corresponding results for AEcap are 
shown in the inset of Fig. [3^ as a function of the incident 
energy. The deviation of the datum point at lowest en- 
ergy from the straight line in the main panel transforms 
into a pronounced increase of AEcap at low energy. 

^From the few times at which we could capture the for- 
mation of a bridge with the camera, we could estimate the 
bridge volume assuming an axisymmetric shape. For the 
bridge at the slowest impact, we obtained V = 0.25±0.06 
for the dimensionless volume. Together with eq. ^ this 
leads to AEcap — 0.214 ± 0.056/j.J, where we have as- 
sumed complete wetting (0 = 0). This result is indicated 
by the shaded area in the inset of Fig. [3]i and compares 
very favorably with the experimental value obtained at 
low impact energy. 

The reduction of AEcap at larger impact energies may 
be understood when one considers the dynamics of for- 
mation of the bridge. There is some time needed for the 
liquid to rearrange from the thick wetting layer around 
the contact point into the liquid capillary bridge struc- 
ture. The entire process spans several time scales, start- 
ing from the microsecond range at individual asperities 
of the grain roughness |28j to ripening processes on the 
scale of minutes [29] . Consequently, although the impact 
duration of the glass sphere with the bottom is sufficient 
to form a small bridge, the volume of the latter will be 
larger when it has more time to form. The characteris- 
tic time scale of the early stage of bridge formation can 
be estimated from the typical height of the bridge at 

contact, which is h = h/R ~ \Jv /2tt « 0.2. Know- 
ing that viscous damping is of minor importance, we 
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FIG. 3: (a) A 'return map' of the total energies between the 
impacts. The accuracy of the measurement is sufficient for 
the deviation of the datum point at smallest energy from the 
line to be significant. It points to finite-time efi'ects in the 
dynamics of bridge formation and rupture, (b) The energy 
loss associated with bridge rupture as a function of bridge 
volume, which was determined independently from images 
taken around the impacts. The solid curve represents the 
quasi-static limit (eq. Q), without fitting parameters. Al- 
though this equation relies on considerable approximations, 
the deviation of the data points at smaller volume is clearly 
significant. 



consider the dispersion relation of undamped capillary 
waves, uj^ = jq^/pnq- With g « 1/h we obtain the time 
scale T — I/lli K 3.3 x 10~^ sec. It is straightforward to 
calculate the kinetic impact energy Er at which the solid 
surfaces are closer than h for a duration equal to r. We 
obtain 



= « 780n J 

ihpuq 



(14) 



which is of the same order as the transition seen in the 
inset of Fig. [3^. 
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However, the dynamics not only of bridge formation, 
but also of bridge rupture gives rise to variations in 
AiScap- If the solid surfaces withdraw rapidly, the for- 
mation of the extended neck will be impeded and the 
bridge is expected to pinch off at a separation which is 
smaller than the 'quasi-static' Sc- This reduces the up- 
per limit of the integral in eq. (|3| , and thus the value of 
Ai?cap- By measuring the bridge volumes for different 
impacts independently, we can distinguish these two ef- 
fects. This is shown in Fig. [3]3, where AEcap is plotted 
as a function of the bridge volume, as determined from 
images close to the respective impacts. The solid line 
represents eq. ([6| and has no fitting parameters. Good 
agreement is found for the larger bridge volume, which 
corresponds to the leftmost point in the inset of Fig. [3^. 
At smaller volume, AEcap is indeed reduced, but this re- 
duction in is much stronger than predicted by the solid 



line, which represents the quasi-equilibrium shapes. 

If we finally return to the conversion of the energy 
loss into an energy-dependent restitution coefficient, we 
see that the variations in AEcap will cause e{E) to de- 
crease even stronger with decreasing impact energy than 
suggested by eq. ([7|. This stresses again the quahtative 
differences to the dry systems, where e{E) tends to be de- 
creasing with increasing impact energy. When one looks 
at the results in the main panel of Fig. [3^, one might not 
anticipate that the small negative intercept with the ver- 
tical axis should be of any importance for the collective 
behavior of many spheres. Quite surprisingly, the dra- 
matic mechanical differences between dry and wet sand 
show that this is nevertheless the case. 
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